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A  NUMERICAL  MODEL  OF  THE  EFFECTS  OF  PROPELLER  WASH  AND 
SHIP-INDUCED  WAVES  FROM  COMMERCIAL  NAVIGATION  IN  AN  EXTENDED 
NAVIGATION  SEASON  ON  EROSION,  SEDIMENTATION,  AND  WATER  QUALITY 
IN  THE  GREAT  LAKES  CONNECTING  CHANNELS  AND  HARBORS 


INTRCDUCriCN 

Ship  traffic  in  river  channels  generates  forces  c jver  and  aboe  the 
natural  forces  created  by  fluid  flew  in  the  channel  that,  when  of  a  sufficient 
magnitude,  can  have  an  adverse  impact  on  the  natural  fluvial  environment. 

These  forces  are  developed,  in  part,  as  a  result  of  propeller  jet  and 
displacement  or  backwater  flow  velocities.  Other  important  forces  are  ship 
waves  and  so-called  "drawdown"  effects.  All  of  these  forces  interfere 
although  some  are  strong  functions  of  channel  section,  vessel  speed,  etc. 
Estimation  of  the  magnitudes  of  these  forces  and  their  resulting  environmental 
effects,  principally  channel  bed  erosion  rates  and  the  increase  of  suspended 
solids,  is  an  important  element  in  the  development  of  plans  for  improvement 
and  optimum  utilization  of  waterways.  The  following  presentation  encompasses 
the  effects  of  propeller  jet  and  backwater  velocities  and  ship- induced  waves. 

Predictive  equations  based  on  a  combination  of  theory  and  observations 
have  been  developed  for  estimating  the  unknown  quantities.  The  equations  are 
used  to  model  forces  created  by  vessels  transiting  harbors  and  interrxrinecting 
channels  of  the  Great  Lakes.  The  mathematical  model  has  been  computerized  in 
FORTRAN  77  code  for  real-time  execution  on  an  IBM  PC-AT  computer  or  a 
compatible  machine  with  at  least  10  megabytes  of  randan  access  memory.  A 
principal  objective  of  the  computer  model  is  to  permit  an  assessment  of  the 
environmental  effects  of  vessel  operations  in  an  extended  navigation  season 
when  ice  cover  is  present.  Example  results  are  presented  for  two  locations  an 
the  St.  Marys  River.  While  the  model  has  been  calibrated  for  the  St.  Marys 
River,  it  is  sufficiently  general  to  be  of  use  in  other  similar  geographic  and 
environmental  settings. 
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ENVIFOWENTAL  SETTING 


The  St.  Marys  River,  an  important  ccirmercial  navigation  artery,  connects 
Lake  Superior  and  Lake  Huron  along  the  border  between  the  United  States 
(Michigan)  ard  Canada  (Ontario).  Water  elevation  drops  about  22  feet  along 
the  60  mile  length  of  the  river.  Much  of  the  elevation  difference  occurs  at 
the  St.  Marys  Falls  which  are  negotiated  by  ships  in  one  step  via  a  set  of 
parallel  locks.  The  river  consists  of  a  succession  of  natural  and  man-made 
channel  sections  which  exhibit  a  variety  of  physical  and  natural 
characteristics  and  dimensions.  In  winter  the  St.  Marys  River  accumulates  an 
extensive  cover  of  ice. 

Several  lines  of  evidence  suggest  that  propwash  from  carmercial  vessel 
passages  causes  a  disturbance  of  St.  Marys  River  channel  sediments.  The 
evidence  includes  observations  of  increased  turbidity  and  suspended-sediment 
concentrations  following  vessel  transits  as  well  as  a  paucity  of  benthic 
organisms  on  the  channel  bed  in  sane  areas,  which  may  not  be  due  to  substrate 
or  temperature  differences.  Disturbance  of  channel  sediments  may  be  enhanced 
during  the  winter  season  when  navigation  in  ice  requires  the  use  of  increased 
vessel  power. 

The  extent  and  intensity  of  disturbance  of  bottom  sediments  is  directly 
related  to  the  nature  of  sediments  present.  Based  on  information  provided  by 
the  Detroit  District,  Corps  of  Engineers,  St.  Marys  River  sediments  range  from 
non-cohesive,  gravel-size  material  to  dense  clays  that  are  strongly  cohesive. 
These  changes  may  be  observed  both  in  the  axial  direction  and  transversely  at 
individual  river  cross  sections. 

During  the  period  mid-December  to  mid-January  ice  begins  to  accumulate  on 
the  St.  Marys  River.  As  the  season  progresses,  the  ice  consolidates  into  a 
solid,  continuous  cover.  At  ary  stage  in  its  development,  the  ice  cover  may 
be  broken  by  vessel  movement.  Conmercial  vessel  activity  normally  results  in 
a  broken  chanrel  through  an  otherwise  continuous  ice  cover.  Icebreaking 
activity  may  yield  a  broken  ice  field  across  the  breadth  of  the  river  channel. 
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MXEL  FORMULATION 


The  numerical  model  to  be  presented  is  based  an  tte  premise  that  in¬ 
channel  farces  resulting  from  vessel  movements  are  created  by  backwater 
velocity,  propeller  jet  velocity  and  ship- induced  waves.  It  is  also  assumed 
that,  to  a  first  approximation,  the  effect  of  water  velocities  induced  by 
vessel  motion  is  of  a  similar  nature  to  the  effect  of  water  velocities 
associated  with  natural  gravity  flews.  The  latter  assumption  provides  the 
basis  for  estimating  relationships  between  sediment  concentration  and  flow 
velocity  and  for  estimating  the  size  and  quantities  of  susficial  bottan 
particles  displaced.  In  this  context,  flew  velocity,  as  used  throughout  this 
report,  is  the  algebraic  sum  of  backwater,  propeller  jet,  and  ambient  river 
velocities. 

General 

As  a  vessel  moves  forward  in  a  restricted  channel,  the  water  displaced  by 
the  movement  flows  in  a  direction  opposite  to  vessel  motion.  Displacement 
flow  velocity  increases  rapidly  as  vessel  speed  approaches  the  so-called  first 


critical  velocity,  formulated  by  Hochstein  (1980)  as  follows 


V  -  KCgA^BJ- 


where 


A  =  the  cross-sectional  area  of  the  channel  (ft) 
c 

B  =  the  surface  width  of  the  channel 
c 

g  =  gravitational  acceleration  (32.2  ft/s/s) 

K  *  the  oonstrairment  factor  defined  as  a  function  of  Va  blocking  ratio 

n  *  A  /A  and  the  ratio  L/b  where  Am  is  the  cross-sectional  area  of  the 
cm  m  m 

submerged  parti  on  of  the  vessel  and  L  and  bm  are  vessel  length  ard  beam, 
respectively.  The  curve  which  defines  K  is  given  by  Hochstein  (1980).  The 
constant  K  <  1.0  in  all  cases. 
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As  V  decreases,  vessel  speeds  must  also  be  decreased  if  adverse  effects 

GT 

are  to  be  minimized  or  prevented.  Equation  (1)  shows  that  decreasing  either  K 
or  the  ratio  A^yB,  (hydraulic  depth)  decreases  Vff.  In  both  cases,  the  result 
is  a  decrease  in  the  value  of  the  blocking  ratio. 

Closely  associated  with  displacement  flew  is  the  phenomenon  referred  to 
as  "drawdown. "  As  the  result  of  vessel  movement  in  restricted  channels,  a 
displacement  wave  of  depression  fonts,  its  phase  velocity  being  equal  to  the 
speed  of  the  ship  which  creates  it.  At  a  ship's  speed  equal  to  the  first 
critical  velocity,  the  displacement  wave  becomes  a  classical  "hydraulic  jump. " 
The  wave  can  be  effective  in  creating  negative  impacts  in  shallow  water  ard  at 
channel  banks. 

Disturbances  also  are  caused  by  the  propeller  jet  which  may  impinge  on 
the  channel  bed  and  banks.  The  jet  is  in  the  form  of  a  cone  in  which  the 
velocitv  rapidly  decreases  awav  from  the  propeller.  Bottom  disturbance,  thus, 
is  a  function  of  the  distance  from  propeller  axis  to  the  bottom.  Propeller 
jet  disturbances  are  greatest  at  bends  in  a  channel  and  in  reaches  where 
cross-winds  may  force  vessels  to  move  with  a  substantial  angle  to  the  channel 
axis. 

When  a  ship  moves  over  a  water  surface,  a  system  of  bow  waves  and  a  wake 
are  formed.  The  bew  waves  ard  wake  are  propagated  in  the  form  of  four  trains 
of  diverging  and  transverse  waves  formed  at  the  bew  and  stem,  respectively. 
The  midlire  of  the  front  of  the  diverging  waves  makes  an  angle  of  about  20 
degrees  with  the  vessel  trackline.  The  transverse  waves  are  included  inside 
the  diverging  wave  ard  their  propagation  (phase)  velocity  is  equal  to  ship’s 
speed  as  is  the  displacement  wave.  In  shallow  water,  orbital  velocities 
associated  with  these  surface  gravity  waves  frequently  are  adequate  to 
resuspend  sedimentary  material. 
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Backwat  Velocity 

Backwater  flew  velocity  is  calculated  using  an  approach  developed  by 
Hochstein  (1967)  which  relates  the  change  in  flew  velocity  A to  vessel  speed 
V  ,  blocking  ratio  and  critical  velocity,  vie. 

iVw  =  Vs[(aB-B>l)°*5-l]  (2) 

and 

-W:V  3  =  maxi1'  °-114(Bc/bm)  +  °-715)  <3> 

and  is  defined  as  the  change  in  mean  strearwise  veiocit/. 


In  Equation  (2),  a  =  (n/(n-l)2‘J  and  B  »  O.Se1 '  ^s^cr,  V  /V  <  0.65  or 

s  cr  — 

B  =  1.,  0.65  <  V  /V  <  1.  (HOchstein,  1967). 


The  shape  of  the  flow  velocity  cross -sectional  diagram  can  be  represented 


by  the  following  function 

vu(y)  >  kx  e('y/V  (4) 

where  v 

V  (y)  =  the  backwater  flow  velocity  (fts-1)  at  distance  y  (ft)  from 


the  sailing  line 


*2 


"  v<°>  "  AV<*>max 
=  Bx/ct  (l-e^F(a)) 


aAv  (i.e.,  the  velocity  at  the  sailing  line) 


B1  =  distance  (ft)  from  vessel  centerline  (sailing  line)  to  shoreline 
F(cc)  =  0.42  +  0.52  In  a 

The  model  yields  an  evaluation  for  any  point  on  the  channel  cross  section.  In 
the  examples  provided  in  the  following  section  it  is  assumed  that  the  vessel 
operates  on  the  chanre  sailing  line  and  never  exceeds  the  established  speed 


limits. 
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The  relative  accuracy  of  the  backwater  equation  far  predicting  backwater 
velocities  is  shewn  in  Figure  1.  The  data  far  Figure  1  were  generated  during 
a  test  of  the  prediction  technique  cn  the  Kanawha  River  in  West  Virginia.  The 
test  was  conducted  under  open-water  ccnditicns.  Conditions  cn  the  Kanawha 
River  where  the  observations  were  made  (channel  configuration,  depth,  etc. ) 
are  not  dissimilar  from  those  of  the  St.  Marys  River. 

Drawdown 

Calculations  of  drawdown  heignts  are  based  cn  the  approach  of  Hochstein 
(1967)  which  uses  Bernoulli's  equation  as  the  starting  point.  With  this 
approach,  the  total  difference  in  elevation  between  the  undisturbed  water 
surface  and  the  lowest  point  cn  the  displacement  wave  is  given  as 

-H  =  Vs2(a-l)B/2g  (5) 

The  relative  accuracy  of  Equation  (5)  for  predicting  drawdown  heights  is 
shown  in  Figure  2.  The  observations  in  Figure  2  were  extracted  from  drawdown 
data  collected  at  seven  cross  sections  along  the  St.  Marys  River,  and 
represent  84  usable  data  points  from  a  total  of  94.  The  data  represent  open 
water  conditions. 

Initial  results  using  Equation  (5)  shows  that  drawdown  values  were 
consistently  underestimated  for  the  Sand  Island  and  West  Cell  Dock  data. 
Analysis  of  all  cross  sections  showed,  far  these  two  sections,  a  relatively 
narrow  channel  section  within  a  much  broader  river  cross  section.  A  dramatic 
improvement  of  the  fit  of  the  data  from  these  two  sections  was  achieved  by 
multiplying  the  calculated  AH  by  the  term  3.4  -  (5.8  x  C_/Bc)  where  Cc  is  the 
width  of  the  maintained  channel  and  C^/B^  £  0.2.  The  data  in  Figure  2  have 
been  so  corrected. 
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Figure  1. — Correlation  plot  of  measured  versus  calculated 

backwater  velocities. 
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Figure  2.  — Correlation  plot  of  observed  versus  calculated  drawdown  heights 
using  84  data  points  from  seven  cross  sections  on  the  St.  Marys  River. 


Sane  cannents  concerning  the  quality  of  the  observational  drawdc*  n  data 
is  essential.  The  level  of  accuracy  at  which  necessary  data  was  recorded 
showed  sane  inconsistencies.  For  example,  in  sane  cases  vessel  draft  was 
recorded  to  the  nearest  0.1  foot  or  even  0.01  foot,  whereas,  in  others  it  was 
only  recorded  to  the  nearest  whole  foot.  Tabulated  speeds  in  at  least  two 
instances  were  different  than  the  values  that  would  be  obtained  by  dividing 
vessel  length  by  the  difference  in  time  for  bow  and  srtem  passage  past  the 
recording  point.  Specifically,  the  speed  of  the  Edgar  B.  Speer  during  a 
dewnbeund  passage  on  5/28/35  was  given  as  9.65  ft/s,  whereas,  the  value 
calculated  using  vessel  length  and  time  difference  was  9.84  ft/s.  For  the 
downbeund  passage  of  the  Canadoc  on  5/29/85,  tabulated  and  calculated  speed 
values  are  12.54  ft/s  and  12.60  ft/s,  respectively. 

Perhaps  the  most  serious  shortcoming  of  the  observational  data  was  the 
failure  to  indicate  ship  position  relative  to  the  various  observation  points. 
An  evaluation  of  the  effect  of  ship  distance  an  wave  height  at  the  channel 
bank,  thus,  was  impossible.  It  is  assumed  here,  therefore,  that  the  vessels 
navigate  always  cn  the  sailing  line.  Moreover,  no  attempt  was  made  during  the 
observation  period  to  separate  the  effects  of  drawdown,  ship  divergent  waves 
and  wind-geivrared  surface  waves.  Calculations  here  and  elsewhere  suggest 
that  these  complicating  factors  may  be  important. 

Given  the  deficiencies  in  the  data,  the  fit  of  observations  versus 
calculations  shewn  in  Figure  2  is  quite  remarkable.  The  goodness  of  fit  is 
shown  statistically  as  well.  The  means  and  standard  deviations  of  the 
observed  and  calculated  drav*Jcwn  heights  were  6.1,  4.9  and  5.9,  4.7, 
respectively. 
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Ice  conditions  inpact  cn  the  drawdown  effect  reflected  through  reduction 
in  channel  cross-secticnal  area  and  is  manifested  by  changes  in  the  values  of 
a  and  B  in  equation  (5). 

The  relative  accuracy  of  equation  (5)  to  predict  dras*3cwn  heights  under 
ice  based  on  only  seven  data  points  is  shewn  in  Figure  3.  Although  the  sanple 
is  too  small  to  draw  any  statistically  significant  conclusions,  far 
completeness  and  to  provide  conservative  (calculated  values  with  tendency  to 
be  higher  than  measured)  evaluation,  the  model  has  been  corrected  to  reflect 
the  difference  between  observed  and  calculated  drawdown  values.  The  least 
squares  line  through  the  origin  shown  cn  Figure  3  has  a  slope  of  1.7  and  this- 
is  the  multiplier  of  equation  (5).  Model  values  for  drawdown  under  ice  should 
be  used  with  caution  until  additional  data  becomes  available  to  verify  the 
validity  of  the  correction  factor. 

Propeller  Jet  Flew  Velocity 

Calculations  of  propeller  jet  flow  velocities  are  based  on  the  approach 
of  Blaauw  and  Van  de  Kaa  (1978)  which  assumes  a  Gaussian  distribution  of 
velocity  within  the  ocne  of  flow  behind  the  propeller. 

The  maximum  propeller  jet  velocity  at  a  radial  distance  r  from  the  axis 
is  given  as 

V  .(r)  -  1.285  S  (D  P  /S)1/3/r  (6) 

pj  p  e 

where 

P  »  vessel  horsepower 

S  »  the  number  of  propellers 

Dp  *  the  propeller  diameter  (ft) 

Figure  4  shows  the  ability  of  the  method  to  predict  propeller  jet  {. 

velocities.  These  data  were  collected  at  the  same  time  and  under  similar 
conditions  to  those  shown  in  Figure  1. 
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OBSERVED  (tn) 


8  10  12 

CRLCUIATED  (in) 


Figure  3.-- Correlation  plot  of  observed  versus  calculated  drawdown  heights 
using  seven  data  points  (two  are  identical)  from  two  cross  sections  on  the 
St.  Marys  River.  Data  collected  with  ice  cover  present.  The  hccken  line 
is  the  least  squares  best  fit  of  the  points. 
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OBSERVED  (NORMALIZED) 


CALCULATED  (NORMALIZED) 


Figure  4. --Correlation  plot  of  measured  versus  calculated 
propeller  jet  velocities. 
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Ice  Operations 


Power  requires  for  vessel  operation  at  specified  speed  limits  is  arrived 
at  with  an  equation  of  the  form 

Pa  =  AV^  (7) 

e  s 

where 

A  and  b  =  constants 

Using  an  example  given  by  Wuebben,  et.al.  (1984)  based  on  a  model  of  the 
MV  ST.  CLAIR,  the  values  for  the  constants  in  equation  ,  ere  found  by  a 
least  squares  technique  to  be  A  =  2.4  and  b  =  3.  The  va_  .  of  course,  are 
applicable  cnly  to  open-water  operations. 

In  ice,  power  required  to  attain  a  given  speed  of  advance  is 
substantially  greater  than  is  given  by  equation  (7).  To  accomodate  this 
difference,  curves  of  vessel  speed  versus  ice  thickness  (Figures  5  and  6) 
provided  by  the  Detroit  District,  Corps  of  Engineers,  were  used  to  scale  the 
constants  in  equation  (7).  The  curves  were  partitioned  into  N  straight  line 
segments,  with  N  sufficiently  large  to  completely  describe  the  behavior  of  the 
curve. 

For  operations  in  ice,  it  was  assumed  that  propeller  jet  velocity,  far 

the  same  expenditure  of  horsepower,  Increased  as  the  square  of  the  velocity 

for  open-water  operations.  The  factor  applied  for  ice  operations  thus  is 

2 

(max.  vessel  speed  in  open  water/vessel  speed  in  ice)  and  serves  to  correct 
both  backwater  and  propeller  wash  velocities.  The  maximum  propeller  jet 
velocity,  of  course,  is  limited  by  the  available  vessel  power  as  shewn  by 
equation  (6).  In  the  example  cases  to  be  presented,  ice  thicknesses  prevented 
vessels  fran  attaining  established  speed  limits  even  with  maxiirum  power 
expenditure. 

Where  ice  cover  is  present,  channel  cross  section  is  effectively  reduced. 

This  reduction  will  result  in  an  increase  in  K,  V  V  ,  &H,  and  V '  as  is 

cc  w  pj 

shown  in  Equations  1  through  6. 
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MnXIMUM  VESSEL  SPEED  (MPIIl 


Figure  5 . ~ Maximal  vessel  speed  of  advance  versus  thicJcvess  of  broken  ice. 


MAXIMUM  VESSEL  SPEEO  (MPM 


Figure  6. — Maxirrun  vessel  speed  of  advance  versus  thickness  of  solid  ice. 
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Under  ice  conditions,  vessel  waves  are  suppressed  and  additional  friction 
is  introduced.  Possible  inpact  of  these  changes  is  complex  and  might  be 
accurately  described  cnly  based  on  intensive  field  measurements.  In  the 
model,  iupact  of  ice-ccnditicns  is  reflected  by  empirical  coefficient  in 
Equation  (5)  for  drawdown  and  by  reduction  in  channel  cross  section. 

Surface  Waves 

For  ship  operations  in  channels  where  the  blocking  ratio  n  >  n ,  only 

diverging  waves  are  important  (Hochstein  aid  Adams,  1985).  Transverse  waves 

are  of  minimal  impact  and  can  be  neglected.  Thus,  for  the  St.  Marys  River 

with  vessel  classes  1-10,  n  >  15  always,  only  diverging  waves  are  treated. 

The  height  of  divergent  waves  as  given  by  Hochstein  (1967)  is 

h  =  0.0448  V  2(d  /L)°*5a  (8) 

s  s 

where 

d  =»  vessel  draft, 
s 

For  most  cases  of  interest,  length  of  divergent  waves  1  -  ICh  (Hochstein, 
1967).  A  well-known  expression  which  relates  maximum  bottom  wave  orbital 
velocity  to  wave  height,  length  and  period  is  (e.g.  Wiegel,  (1964)) 

Umax  *  a  s/sinh<k  D> 

where 

a  =  wave  amplitude 
s  *  radial  wave  frequency  (2tt/T) 
k  *  wave  number  (2V1) 

D  ■  water  depth 

Because  wave  motions  are  orbital,  wave-associated  velocities  do  not  contribute 
to  net  channel  flow. 
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Sediment  Transoort 


Two  types  of  soil  or  sediment  disturbance  may  te  expected  as  the  result 
of  ship  traffic  -in  restricted  channels.  The  first  type,  referred  to  as  bed 
load,  occurs  when  sedimentary  particles  move  in  proximity  to  the  bottom.  For 
sediment  transport  to  be  classified  as  bed  load,  the  particles  must  be  in 
contact  with  the  bottom,  at  least  part  of  the  time.  Sedimentary  particles 
carried  wholly  in  suspension  and  with  no  bottom  contact  comprise  tlie  suspended 
load. 

Calculations  of  bed  load  transport  are  based  cn  the  method  of  Yalin 
( 1963 ) .  The  Yalin  equation  written  in  terms  of  the  volume  flux  of  sediment 
per  unit  flew  width  is 

Q  =  ax  U*S0[1  -  In  (l+ajS)]  (9) 

where 

.  a1  =  0.635 

32  =  2,45  (c/°s)°-4[vc/(ps-p)g0]0'5 

In  equation  (9),  U*  (=(vq/p )^'5)  is  the  local  shear  velocity,  is 

critical  shear  stress  (to  be  subsequently  discussed)  P  and  P  are  fluid  and 

s 

sediment  densities,  respectively,  and  0  is  the  nominal  grain  diameter.  The 

parameter  S  =  (v  A>  )  -  1  is  the  normalized  excess  shear  stress, 
o  c 

Calculations  of  suspended  load  transport  follow  the  technique  elaborated 
upon  by  Smith  (1977)  and  Adams  and  Weatherly  (1981).  This  method  treats  a 
two-phase  flew  (sediment  and  water)  in  which  the  differential  equation  of  mass 
conservation  represents  a  balance  between  upward  turbulent  diffusion  and 
downward  gravitational  settling.  Integration  of  this  equation  for  a  channel 
flew  situation  yields  an  expresicn  for  sediment  concentration,  viz. 

C  100  6.24p( z-0.2h) 

C  =  xi§-(— )P  exp( - ^ - )  (10) 

o 
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where 

h  ■  water  depth 

p  -  w3/0.4U*z  (B  -  6.24;  Smith  1977) 

z  *  height  above  bottom 

w  *  particle  settling  velocity 

C  =  a  reference  concentration 
o 

The  expression  for  the  reference  concentration  as  developed  by  Smith  and 
McLean  (1977)  is 


where 

i  =  the  fraction  of  the  bed  sample  in  the  median  size  class 

Y  =  2.4  x  10'3. 
o 

As  fluid  velocity  or  bottom  stress  increases  above  a  sediment  bed,  a 
situation  is  reached  such  that  an  additional  increase  causes  sedimentary 
particles  to  be  put'  in  motion.  When  motion  is  incipient,  the  "threshold" 
condition  has  been  met.  Calculations  of  sediment  threshold  are  based  cn  the 
work  of  Miller,  et.al.  (1977)  for  ncn-cohesive  sediments  (gravel,  sand,  coarse 
silt)  and  of  Migniot  (1968)  for  cohesive  materials.  These  expressions  are 

v  =  0.32  Y°*5  0  <  0.08  on  (12) 

c 

where 

Y  ■  shear  strength 

v  -  46.6  00*54  0.08  an  <  0  <  0.2  am  (13) 

c 

v  -  79.4  0°* 9  0  >  0.2  cm  (14) 

c 

It  is  Important  to  note  that  sediments  smaller  than  about  0.02  an 
diameter  will  always  move  in  suspension  as  was  demonstrated  by  Bagnold  (1954) 
and  elaborated  upon  by  McCave  (1971). 
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Model  Parameters  and  Incut 


Certain,  well  defined,  variables  and  parameters  control  the  level  of 
fluid  velocity  induced  by  propeller  jet  and  backwater  effects.  The  values 
assigned  to  each  are  functions  both  of  vessel  (operational  parameters)  and 
channel  (physical  parameters)  characteristics.  The  exact  nature  of  these 
factors  and  the  manner  in  which  they  and  their  variations  are  treated  in  the 
model  have  been  discussed. 

The  Great  Lakes  fleet  vessel  characteristics  required  for  modeling 
purposes  and  the  values  used  in  example  computations  were  provided  by  the 
Detroit  District  and  are  shewn  in  Table  1.  It  should  be  noted  that 
operational  data  for  vessel  operation  in  ice,  specifically  the  speed- 
horsepower  relationship  could  be  obtained  only  for  three  vessel  types,  i.e. 
classes  5,  8,  and  10.  Values  appropriate  to  classes  6,  7,  and  9  were  obtained 
by  interpolating  the  available  data.  Vessel  classes  1-4  were  assumed  to  be 
unimportant  because  they  are  likely  to  be  insufficiently  powered,  are 
uncorman,  and  as  a  result  they  are  not  likely  to  have  much  of  an  environmental 
effect.  Speed-horsepower  relationships  were  available  for  broken  ice 
navigation  as  well  as  for  the  situation  in  which  a  continuous  ice  cover  is 
present. 

The  model  results  focus  cn  two  specific  examples,  Frechette  and  Ninemile 
Points  (herein  after  referred  to  as  Sections  1  and  2).  The  appropriate  values 
of  physical  ard  sedimentary  characteristics  of  these  areas  that  are  required 
inputs  to  the  model  are  shown  in  Figures  7  and  8  and  listed  in  Tables  2  and  3. 
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Table  1 . --Characteristics  of  Great  Lakes  Vessels1 


(Model  Input  Values) 


Class 

Length 

Beam 

Propellers 

Number  Diameter 
( feet ) 

Horsepower 

Draft 

( feet ) 

in  feet 

Upbound 

Downbound 

1 

300 

49 

1 

10.5 

3,000 

14 

20 

2 

430 

60 

1 

14.0 

6,000 

20 

22 

3 

524 

60 

1 

14.0 

2,  500 

17 

25 

4 

579 

60 

1 

14.75 

4,000 

15 

26 

5 

620 

68 

1 

17.5 

8,000 

18 

26 

6 

680 

75 

1 

18.5 

7,000 

19 

26 

7 

730 

75 

1 

17.5 

9,600 

20 

26 

8 

767 

75 

1 

18.5 

8,500 

19.5 

26 

9 

9C0 

105 

1 

20.0 

13,600 

20 

26 

10 

1,000 

105 

2 

17.5 

19,500 

24 

27 

1 


Data  provided  by  Detroit  District, 


Corps  of  Engineers. 


Table  2. — Physical  characteristics  of  St.  Marys  River  cross  section 


Location 

Frechette  Point 

Ninemile  Point 

(Section  1) 

(Section  2) 

Upbcund - 

-  11.7 

14.7 

Speed  limit  ( ft/s ) 

Dcwnbound - 

-  14.7 

14.7 

Ambient  velocity  (ft/s) - 

-  2.7 

0.8 

Ice  thickness  (ft) - 

-  0.88 

1.11 

Table  3. — Characteristics  of  St.  Marys  River  sediments 


Frechette  Point  Ninemile  Point 

Location - - 

(Section  1)  (Section  2) 


Position 

Description 

Diameter  (mm) 

Description 

Diameter  (im) 

Left  bank - 

Clay 

0.004 

Clay 

0.002 

Mid-channel — 

Clay 

0.004 

Sand 

0.300 

Right  bank — 

Clay 

0.004 

Silt 

0.005 

20 


DEPTH  (FT) 


FRECHETTE  SECTION  (1) 


NINE  Mile  POINT  SECTION  (2) 


0  10  20  30  UO  50  60 

WIDTH  (FT)  «102 


Figure  8. — Ninemile  Point  cross  section  (Section  2). 
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MODEL  RESULTS 


The  results  presented  below  are  based  on  computations  using  the  model 
described  above.  Specific  outputs  of  the  model  include  propeller  jer, 
backwater  and  net  velocities,  suspended-sediment  concentration,  bed  load 
sediment  flux,  kinetic  energy  density,  drawdown  and  di'/ergent  wave  height. 
Kinetic  energy  density  is  here  defined  as  cne-half  times  the  net-velocity 
squared.  The  first  six  parameters  have  been  calculated  at  eleven  <1"  )  poir.ts 
across  a  channel  section,  at  the  sailing  line  and  at  5  equidistant  points 
either  side  of  the  sailing  line.  Because  the  navigation  channel  frequently  is 
nearer  one  bank  than  the  other,  partitioning  the  cress  section  in  this  manner 
as  opposed  to  equal  spacing  of  points  across  the  entire  channel  width 
appropriately  weights  impacts  to  that  side  of  the  channel  cm  which  they  are 
created.  Sailing- line  values,  calculated  for  up  bound  and  dewnbound  navigation 
fer  all  10  vessel  classes  in  open  water  and  6  vessel  classes  in  ice  are  shewn 
in  Tables  4  and  5  for  Sections  1  and  2,  respectively.  Because  of  the  lew 
power  assigned  to  vessel  classes  1-4  (Table  1),  it  has  been  assumed  that  they 
will  be  unable  to  operate  in  ice.  For  completeness,  however,  computations 
with  respect  to  vessel  classes  1-4  have  been  made  and  the  results  are  shewn  in 
all  subsequent  data  presentations.  The  model  results  indicate  that  vessel 
impacts  with  respect  to  kinetic  energy  density,  suspended-sediment 
concentration  and  sediment  flux  are  concentrated  at  the  sailing  line.  The 
results  to  be  presented  subsequently  are  for  the  sailing  line  unless  otherwise 
indicated.  Furthermore,  suspended- sediment  concentrations  represents  an 
average  of  three  values,  near-bottom,  near-surface,  and  mid-depth. 
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It  is  Important  to  note  that  the  accuracy  of  model  output  values  are 
strongly  dependent  on  the  quality  of  input  parameters,  particularly  in  the 
case  of  turbidity  are!  suspended  sediments.  These  outputs  must,  therefore,  be 
viewed  in  the  context  of  parametric  inputs. 

Shown  in  Figure  9  is  the  kinetic  energy  density  associated  with  each 
vessel  class  for  open-water  and  ice  navigation  through  Section  1.  In  this  and 
subsequent  portrayals,  ICE  1  refers  to  broken  ice  while  ICE  2  represents  a 
solid  ice  sheet.  The  values  shown  are  for  the  sailing  line.  Kinetic  energy 
density  decays  to  not  more  than  6  percent  of  the  sailing  line  value  at  the 
first  grid  point  an  either  side  of  the  sailing  line.  Thus,  almost  95  percent 
of  the  kinetic  energy  associated  with  ship  passage  is  confined  to  the  center 
of  the  channel. 

Several  interesting  aspects  of  the  problem  emerge  from  Figure  9.  In  each 

case,  i.e.  open  water,  ICE  1  or  ICE  2,  the  kinetic  energy  per  vessel  passage 

is  greater  for  an  upbound  (enpty)  transit  than  for  a  dcwnbcund  (loaded)  one. 

Also,  kinetic  energy  densities  associated  with  navigation  of  vessel  classes 

5-10  in  ICE  1  (broken)  are  greater  than  those  associated  with  operations  in 

ICE  2  (sheet)  ice.  Finally  and  perhaps  most  significantly,  the  magnitude  of 

the  energy  effects  of  vessel  class  10  is  substantially  greater  than  the  other 

nine  classes.  For  example,  kinetic  energy  associated  with  an  open  water, 

upbound  passage  of  a  class  10  vessel  is  greater  than  that  for  any  other  vessel 

operating  in  ice,  regardless  of  the  direction  of  travel.  The  magnitude  of 

kinetic  energy  effects  of  a  class  10  vessel  are,  in  the  worst  case,  about 

* 

twice  as  great  as  those  of  any  other  vessel  class.  The  greater  forces 
generated  by  vessel  class  10  are  apparently  explained  by  at  least  three 
reasons:  (1)  the  only  vessel  type  with  2  propellers,  (2)  high  horsepower  and 
ability  to  maintain  high  speed  in  ice  conditions,  and  (3)  a  draft  of  27  ft. 
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Figure  9. — Curves  of  kinetic  energy  (density  versus  vessel  type.  Section 

at  the  sailing  line. 


27 


Examination  of  Figure  10  which  is  kinetic  energy  density  at  the  sailing 


lira  associated  with  navigation  through  Section  2  yields  results  somewhat 
different  from  those  found  for  Section  1.  For  example,  the  magnitudes  of 
kinetic  energy  are  lower  for  most  vessel  classes.  This  can  be  directly 
related  to  the  greater  cross-sectional  area  of  Section  2.  For  vessel  classes 
1,  2,  7,  8,  and  10,  in  open  water  kinetic  energy  per  vessel  passage  is 
greatest  for  an  upbound  vessel,  whereas,  for  vessel  classes  3-6,  and  9,  it  is 
greatest  for  a  dcwnbound  vessel.  For  a  class  10  vessel,  kinetic  energy 
effects  for  upbound  and  dcwnbound  vessels  in  open  water  are  the  same.  In  ice, 
kinetic  energy  per  vessel  passage  for  all  classes  is  greater  for  a  downbcund 
transit.  While  the  absolute  magnitudes  at  Section  2  are  smaller  than  at 
Section  1  the  relative  effect  of  a  class  10  vessel  is  somewhat  greater  at 
Section  2  (i.e.  the  worst  case  situation  at  Section  2  yields  a  ratio  of 
approximately  3). 

As  noted  above,  at  Section  1,  kinetic  energy  per  vessel  passage  is  always 
greater  for  an  upbound  then  for  a  dcwnbound  vessel.  At  Section  2,  however, 
certain  vessel  classes  induce  a  greater  affect  when  dcwnbound  in  open  water; 
dcwnbound  passages  in  ice  induce  a  greater  affect  than  upbound  passages  focr 
all  vessel  classes.  The  reason  for  these  differences  is  due  to  the  difference 
in  cross-sectional  areas  of  the  two  sections  and  the  associated  ambient 
velocities.  For  upbound  vessels,  propeller  wash,  backwater,  and  ambient 
velocities  have  the  same  sign,  whereas,  for  downbcund  vessels,  propeller  wash. 
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and  backwater  velocities  are  of  opposite  sign  to  that  of  ambient  velocity. 
Backwater  and  propeller  wash  velocity  are  greater  for  downbound  vessels  as  the 
blocking  ratio  is  snail er  and  the  propeller  axis  is  nearer  the  channel  bed 
because  of  the  greater  draft.  Thus,  at  Section  1,  the  effects  of  upbound 
passages  of  class  5  vessels  will  always  be  greater  than  dcwnbcund  passages  as 
long  as  ambient  velocity  is  greater  than  1.5  ft/s.  At  Section  2,  because 
ambient  velocity  is  less  than  1.0  ft/s,  a  downbound  passage  of  a  class  5 
vessel  will  induce  a  greater  affect  than  an  upbound  passage.  Similar 
reasoning  applies  to  the  other  vessel  classes  as  well. 

The  absolute  Impact  of  downbound  (loaded)  vessels  is,  however,  always 
higher  due  to  associated  smaller  blocking  ratio  and  similar  vessel  speeds. 

This  is  manifested  in  higher  drawdown  values  far  downbound  vessels. 

Suspended-sediment  concentrations  induced  at  the  sailing  line  at  Section 
1  by  the  passage  of  the  various  vessel  classes  are  shown  in  Figure  11.  The 
relationships  between  suspended-sediment  concentration  and  vessel  class  are 
similar  to  those  between  kinetic  energy  density  and  vessel  class  shewn  in 
Figure  10.  Suspended-sediment  concentrations  induced  by  the  passage  in  open 
water  of  vessel  classes  1-9  range  from  0-7  mg/L.  When  ice  is  present  these 
values  may  increase  slightly.  A  maximum  of  10  mg/L  can  be  associated  with 
vessel  class  9  in  broken  ice.  The  effect  of  vessel  class  10  is  to  create  a 
suspended- sediment  concentration  2-3  times  greater  than  that  resuspended  by 
other  vessel  types.  The  difference  in  sediment  concentration  between  open- 
water  and  ice  navigation  for  class  10  is  approximately  20  mg/L. 
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Figure  11. — Curves  of  suspended- sediment  concentration  versus  vessel  type. 

Section  1,  at  the  sailing  line. 
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At  section  2,  the  bottom  remains  undisturbed  far  all  vessel  classes  in 
open  water  and  in  sheet  ice.  Passage  of  a  class  10  vessel  through  broken  ice 
resuspends  cnly  trace  quantities  of  sediment  (Table  5).  The  difference  in 
sediment  resuspensicn  at  Sections  1  and  2  is  due  not  cnly  to  the  difference  in 
cross-secrticnal  areas  tut  also  to  the  prevalence  of  cohesive  sediments  at 
Section  1.  It  should  be  recalled  that  the  suspended- sediment  concentrations 
shown  in  Figure  11  are  depth  weighted  averages  of  three  values  (near  bottom, 
near  surface,  and  mid-depth). 

Passages  of  class  1-10  vessels  through  Section  1  failed  to  create  any  bed 
load  transport.  This  result  is  expected  because  of  the  absence  of  particles 
of  a  size  that  would  normally  move  as  bed  load.  Although  the  preserve  of  sand 
at  Section  2  makes  bed  load  transport  possible  the  model  predicts  only  trace 
amounts  of  sediment  flux  induced  by  class  10  vessels  operating  in  broken  ice 
(Table  5).  For  upbound  transits,  sediment  would  be  displaced  in  a  downstream 
direction.  For  dcwnbound  transits  sediment  moves  upstream. 

Drawdown  heights,  as  calculated  at  the  bank,  are  shewn  in  Figures  12  and 
13  for  Sections  1  and  2,  respectively.  The  curves  in  both  figures  are  similar 
in  that  vessel  classes  5-8  create  drawdown  waves  of  similar  height,  the  waves 
created  during  ice  operations  being  about  90  percent  higher  than  those 
generated  during  open  water  passages.  The  magnitude  of  drawdown  in  broken  and 
solid  ice  is  about  the  same.  The  maximum  height  of  the  drawdown  waves  for 
vessel  classes  5-8  is  about  0.6  foot  in  Section  1  but  cnly  about  0.3  foot  in 
Section  2.  Drawdown  elevatiens  are  abcut  twice  as  large  for  vessel  classes  9 
and  10  than  they  are  for  the  others.  Here  again  the  values  associated  with 
ice  passages  are  about  90  percent  greater  than  those  for  open  water 
operations . 
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Figure  12. — Curves  of  drawdown  versus  vessel  type.  Section  1. 
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Figure  13. --Curves  of  drawdcwn  versus  vessel  type.  Section  2. 
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A  significant  difference  between  this  model  ard  an  earlier  version  is  the 
incorporation  of  the  effects  of  diverging  waves  in  the  current  one.  A  result 
of  this  additional  feature  is  the  ability  to  estimate  resuspensicn  of 
sediments  at  the  banks  by  orbital  motions.  Model  results  stow  that  divergent 
waves  do  indeed  resuspend  sediments  in  the  shallower  parts  of  the  cross 
section.  In  the  cases  tested  here,  resuspended  materials  were  never  mere  than 
2  mg/'L  and  frequently  were  substantially  less.  It  is  difficult,  however,  to 
i  sold  in  the  effects  of  divergent  waves  from  those  of  drawdown.  It  is  possible 
that  many  of  the  impacts  attributed  to  drawdown  may  be  due  to  divergent  waves 
or  perhaps  to  a  ccmbination  of  the  two. 

To  summarize  the  effects  of  vessel  passages  on  the  resuspensicn  of 
channel  sediments  it  is  sufficient  to  consider  vessel  classes  5  ard  10  only  at 
Section  1.  Shewn  in  Table  6  are  the  relative  increases  of  suspended  sediment 
for  class  5  and  10  vessels  as  a  result  of  the  presence  of  the  two  different 
ice  types.  The  downbound,  ice  free  ccrditicn  is  taken  as  unity. 

Examination  of  the  data  of  Table  6  yields  sane  interesting  results  with 
regard  to  the  different  influences  of  the  two  different  ship  types  cn  the 
cross  section  studied.  Perhaps  the  most  significant  point  is  that  the  range 
of  relative  increase  of  suspended  sediments  for  the  three  variables 
considered,  i.e.  vessel  type,  vessel  direction,  and  ice  type  is  substantially 
different,  varying  from  0  to  25,  the  latter  value  associated  with  an  upbound 
class  10  vessel  in  broken  ice.  The  smallest  impact  of  a  class  10  vessel  is 
similar  to  the  largest  impact  of  a  class  5  vessel.  In  general,  due  to  the 
small  Increase  in  water  velocity  generated  by  a  vessel  relative  to  ambient 
water  velocity,  upbound  values  are  higher  than  dewnbound  values. 
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Table  6. 

--Relative 

abundances  of  suspended 

sediment 

Section  1 

Ship 

Type 

Ice 

Broken 

Solid 

Free 

Ice 

Ice 

5  Upbound - 

4 

6 

3 

Downbound- 

1 

3 

0 

10  Upbound - 

15 

25 

22 

Downbound- 

7 

17 

14 

To  examine  the  cumulative  effects  of  vessel  traffic  far  an  entire  year, 
the  sum  of  the  product  of  vessel  kinetic  energy  density  times  the  number  of 
vessel  transits  (Table  7)  was  calculated.  The  data  are  plotted  in  Figures  14 
and  15  with  curves  for  the  season  ending  December  15,  for  an  extension  of  the 
season  to  January  1,  and  for  an  extension  of  the  season  to  February  15.  The 
curves  merge  and  are  essentially  indistinguishable.  The  result  of  the  season 
extension,  thus,  is  a  negligible  increase  of  total  energy  density  at  Sections 
1  and  2.  It  should  be  noted  that  the  period  of  December  15-February  15  was 
assumed  to  be  characterized  by  a  broken  rather  than  a  solid  ice  cover.  Higher 
impacts  are  associated  with  the  broken  cover  than  with  the  solid  oover. 


Table  7 .  -  -Vessel  Traffic  Scenario  Extended  Season  Navigation 


1  Apr  -  15  Dec  1979 - 

Vessel  class  breakdcwn 


ten - 

eight - 

seven - 

six - 

five - 

four - 

three - 

two - 

one - 

Direction  to  15  Dec 


upbound - 

downbound - 

16  Dec  -  1  Jan  1980 - 

Vessel  class  breakdown 


ten - 

nine - 

eight - 

seven - 

six - 

five - 

*  two - 

*  one - 

2  Jan  -  15  Feb  1980' 


Passages  through  the  Sault  Locks 


U.S.  Prejected  Canadian  Total 


452 

3,814 

10,266 

586 

1 

2,586 

374 

1 

1,375 

808 

2,287 

3,095 

485 

286 

771 

101 

572 

2,673 

215 

15 

230 

162 

162 

107 

126 

233 

242 

286 

528 

121 

1.767 

5,888 

959 

1,807 

5,766 

100 

43 

143 

12 

12 

4 

4 

18 

1 

19 

7 

2 

9 

10 

1 

11 

47 

2 

49 

22 

22 

2 

15 

/ 

17 

64 

28 

92 

Vessel  class  breakdcwn 


ten -  8  - 

nine -  2  - 

eight -  19  - 

seven -  4  1 

six -  7  - 

five -  24  2 

*  two -  -  14 

*  one -  -  11 


8 

2 

19 

5 

7 

26 

14 

11 


Direction  after  15  Dec 


upbcund -  7 

downbound -  157 


27  34 

44  201 


Vessel  classes  one  and  two  are  included  in  class  five 
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Figure  14. — Total  kinetic  energy  density  versus  vessel  type  for  a  normal 
season  (lower  curve)  for  an  extension  to  January  1  (middle  curve), 
and  for  an  extension  to  February  15,  Section  1. 
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Figure  15. — Total  kinetic  energy  density  versus  vessel  type  for  a  normal 
season  (lower  curve)  for  an  extension  to  January  1  (upper  curve), 
and  for  an  extension  to  February  15,  Section  2. 


Finally,  it  is  of  practical  interest  to  consider  the  duration  over  which 
a  sediment  resuspensian  event  initiated  by  a  vessel  passage  has  a  detectable 
impact  upon  suspended-sediment  concentration  at  a  given  cross  section. 
Hochstein  and  Adams  (1985)  present  a  method  for  estimating  duration  of  impact 
of  resuspension  that  incorporates  vessel  motion  and  particle  settling 
velocity.  The  appropriate  expression  is 

DI  =  L/(  V  +  VJ  +  D(d)  (15) 

S  a 

where 

DI  =  impact  duration 
V  =  ambient  river  velocity 
D  =  duration  of  particle  settling 
d  *  particle  diameter 

and  the  other  terms  are  as  given  previously.  The  particle  settling  velocity, 
of  course,  is  difficult  to  define  in  a  heterogeneous  size  assemblage  of 
particles.  When  the  size  distribution  is  poorly  known  this  problem  becomes 
more  difficult  and  it  is,  perhaps,  better  to  frame  the  problem  in  terms  of 
maxima  and  minima. 

In  the  case  of  Section  1,  data  provided  by  the  Detroit  District,  Corps  of 
Engineers,  indicates  the  bottcm  near  mid-channel  to  be  composed  of  firm  clay. 
CLay  particle  "diameters"  generally  range  from  about  1  to  4  microns. 
Information  provided  by  Hodek,  et.al.  (1985,  p.  AA7)  indicates  that  sediments 
from  near  mid-channel  may  be  as  coarse  as  75  microns.  Settling  velocity  w,  of 
these  particle  sizes  in  water  of  18°C  is  1.36  x  10  cm/s  for  4  micron 
particles  and  4.41  x  10_1  an/s  for  75  micron  particles  (Gibbs,  et.al.,  1971). 
For  these  values  of  w,  the  second  term  in  equation  (15)  above,  thus,  ranges 
from  about  100  hours  to  only  18  minutes  where  water  depth  is  taken  as  16  ft 
(one-half  the  depth  at  Section  1).  The  first  term  is  much  less  than  the 
second  and  for  the  purpose  of  estimating  inpact  duration,  it  may  be  neglected. 


40 


For  very  fine  particles  with  low  settling  velocities,  the  potential  far 
significant  increases  of  suspended- sediment  concentration  is  great.  The 
ahsence  of  a  noticeable  degradation  of  water  quality  during  the  open  water 
navigation  season  suggests  that  if  sediments  indeed  are  eroded  by  the  passage 
of  ships  through  Section  1,  those  sediments  must  be  substantially  coarser  than 
4  microns.  Calculations  based  on  data  from  Table  7,  indicate  that  there  is  a 
ship  passage  about  every  34  minutes  during  the  open  water  navigation  season. 
With  an  impact  duration  of  only  18  minutes,  all  material  resuspended  by  one 
vessel  has  settled  to  the  bottom  before  another  vessel  passes.  For  extensions 
beyond  December  15,  the  difference  between  impact  duration  and  vessel  passage 
period  is  even  greater. 

There  are,  of  course,  no  impacts  at  Section  2  as  materials  are  neither 
resuspended  or  moved  as  bed  load. 

MODEL  CALIBRATION 

The  model  detailed  above  has  been  calibrated  with  the  best  information 
available.  Selection  of  parameter  values  was  made  by  an  analysis  of  data 
available  from  the  St.  Marys  River  and  other  similar  aquatic  environments, 
most  notably  the  Kanawha  River  in  West  Virginia. 

The  basic  equations  in  the  model  have  shewn  good  fit  with  empirical  data. 
Wherever  necessary,  sane  of  the  equations  have  been  modified  or  new  enpirical 
coefficients  introduced,  to  reflect  specifics  of  the  St.  Marys  River 
hydrographic  features.  As  the  result,  the  presented  model  provides  credible 
accuracy  and  sufficiently  describes  magnitude  and  tendencies  of  the  very 
complex  effects  generated  by  a  vessel  movement  in  a  restricted  channel .  The 
rmdel  provides  functional  relationships  between  a  variety  of  variables, 
determining  extent  of  physical  impacts  of  commercial  navigation  during  both 
formal  and  extended  navigation  seasons.  Calibration  has  shown  that  the  model 
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results  by  level  of  accuracy  and  scope,  as  they  relate  to  the  magnitude  of  the 
effects,  are  sufficient  to  make  reliable  conclusions  regarding  extent  of 
rnpacts  due  uo  different  options  of  catineioial  navigation  in  the  Great  Lakes 
connecting  charnels  and  harbors. 

For  the  model  tuning,  however,  it  is  desirable  that  field  data 
corresponding  to  model  output  values  be  collected  in  the  St.  Marys  River. 

These  data  include  but  are  not  necessarily  limited  to  water  velocities 
associated  with  ship  passages  in  open  water  and  under  a  modest  variety  of  ice 
conditions. 

Velocity  data  collection  ideally  would  entail  a  current  measurement 
program  using  self-recording  current  meters  at  one  or  two  points  cn  a  given 
cross  section.  These  devices  obviate  the  need  for  making  measurements  an  the 
ice  under  difficult  and  often  hazardous  conditions.  A  water  sampling  device 
affixed  to  the  above  mentioned  current  meter  would  provide  a  sample  necessary 
to  determine  suspended-sediment  concentration.  Such  an  arrangement  has  been 
successfully  employed  in  sediment  transport  studies  near  the  Mississippi  and 
Atchafalaya  Deltas  cn  the  Louisiana  coast.  Sampling  of  the  bed  load  would  be 
necessary  to  determine  sediment  flux. 

The  sediment  transport  ocmpcnent  of  the  model  requires  a  knowledge  of  a 
number  of  sediment  properties,  all  of  which  require  a  sediment  collection 
program  and  laboratory  analysis  for  their  determination.  The  properties  of 
interest  include  size  distribution,  settling  velocity,  and  yield  strength  (far 
cohesive  sediments).  While  sane  laboratory  analyses  are  essential,  a 
concurrent  parameter  or  sensitivity  study  utilizing  the  numerical  model 
significantly  reduce  the  amount  and  nature  of  the  laboratory  work  required. 

A  sensitivity  analysis  would  be  a  valuable  adjunct  to  the  modeling  work. 

A  properly  formulated  test  would  obviate  some  of  the  work  that  would  normally 
be  undertaken  as  part  of  a  calibration  activity. 
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